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The recent measurement of the muon anomalous mag¬ 
netic moment of the muon by the E821 experiment in 
Brookhaven [Q, 

= 116 592 020(160) x 10"“ (1) 

signihcantly deviates from the Standard Model. Theo¬ 
retical predictions, which include hve-loop QED and two- 
loop electroweak effects, rely on experimental data and 
theoretical models to account for the non-perturbative 
hadronic contributions. Depending on the treatment of 
the latter, the discrepancy between the experiment and 
the theory can be as large as 2.6 cr. After the release of 
the E821 result (p, the hadronic effects came under re¬ 
newed scrutiny and the reliability of various estimates has 
been disputed il- The main focus of those discussions 
was the hadronic vacuum polarization which modifies the 
photon propagator and has been evaluated using data on 
e“'"e“ annihilation into hadrons and the r lepton hadronic 
decays I, § |. 


(a) (b) 

FIG. 1: Pion pole contributions to the muon g — 2. 

Very recently it has been pointed out 0, || that a 
signihcant part of the discrepancy between the theory 


and the experiment may be due to the theoretical eval¬ 
uation of the pion pole contribution to hadronic light- 
by-light scattering, which influences g^ — 2 via diagrams 
shown in Fig. |^. Namely, it has been claimed that al¬ 
though the magnitude of those contributions computed 
in |, |I^, 0, 0 is correct, they had been taken with 

an incorrect (negative) sign. 

Since the pion pole contribution is the largest among 
the hadronic light-by-light scattering contributions to a^, 
the change in sign has important implications for the 
comparison of the experimental result with the theoreti¬ 
cal prediction. If the correct sign is positive, the theory 
and experiment are in much better agreement (about 1 a 
is removed from the reported discrepancy). Motivated 
by this, we have recalculated the pion pole contribution 
to the light-by-light scattering correction to a^. Using 
the vector dominance model for the pion transition form 
factor, we obtain 

^LBL.7r“ ^ ^ (2) 

thereby confirming the result of Ref. |^]. The purpose 
of this Letter is to describe the calculation leading to 
Eq. (^. We do not address the validity of the model. 
Our objective is to obtain an analytical (rather than nu¬ 
merical) result and settle the issue of its sign. 

At low energies, the interaction of a neutral pion with 
photons is described by the Wess-Zumino-Witten La- 
grangian, 

-Cwzw = - p (3) 

where = 3 is the number of colors and ~ 92.4 
MeV is the pion decay constant. Since this is a non- 
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renormalizable interaction, employing it in loop calcula¬ 
tions results in ultraviolet divergences. While this is not 
a problem in principle, since divergent contributions can 
be absorbed into higher dimensional operators of the chi¬ 
ral Lagrangian, in practice this precludes any numerical 
estimate of the corresponding contributions because, as 
in the case of muon anomalous magnetic moment, the 
relevant counterterms are not known. In this situation 
one resorts to models to obtain a finite result. A sim¬ 
ple and commonly adopted option is to introduce a form 
factor into the 71^77 interaction vertex, which damps the 
contributions of highly virtual photons. This results in 
the following 77^77 interaction vertex: 
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where < 71^2 denote the momenta of the two outgoing pho¬ 
tons. 

The transition form factor 77077 ( 91 ) 92 ) depends on 
the adopted model. For our purpose it is sufficient to use 
the simplest of these models (VMD) where it is assumed 
that the transition form factor reads 


, Q 2 ) 


M 2 
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where the parameter M is phenomenologically deter¬ 
mined to be close to the mass of the p meson M « nip ~ 
769 MeV. 

The simple structure of the form factor opens a way 
for analytic calculations. Initially, the diagram has three 
widely separated scales: 


ml - ml ml <s: M^ (6) 


and one can take advantage of this hierarchy to construct 
an expansion in 5 = {ml — ml)/ml and ml/M'^ using 
asymptotic expansions. The expansion in <5 is especially 
simple and reduces to the Taylor expansion of the pion 
propagator. The expansion in m^/M^ is the so-called 
Large Mass Expansion I 100 - 

Consider first the diagram in Fig. |^(a). The lines de¬ 
noted by 7 -propagators are in fact products of the mass¬ 
less photon propagators 1 /kf and of the form factor terms 
fvF /{kf — M^). The vector meson mass M in the latter 
sets the hard scale for the momentum integrals, whereas 
the muon mass sets the soft scale. (After the Taylor ex¬ 
pansion in 6 these are the only mass scales present.) We 
can now obtain four possible combinations of soft/hard 
integration momenta in the two loops of Fig. 0 (a). The 
leading quadratic logarithm of the muon and p mass ratio 
is obtained by adding contributions with both momenta 
soft, both momenta hard, and the one with the momen¬ 
tum in the triangle (left-hand side) loop hard, and the 
right-hand side loop soft. The fourth, soft-hard combi¬ 
nation contributes only in the next order, suppressed by 
m 2 /M 2 . 

The non-planar diagram in Fig. 0(b) can be evaluated 
in a similar manner, with the only difference that a fifth 


momentum region appears. Namely, the two integration 
momenta may be hard, but their difference in the pion 
propagator may be soft. In this case the diagram factor¬ 
izes into a product of two one-loop vacuum diagrams. 

Integrations in all four regions can be carried out ana¬ 
lytically, and coefficients in both expansion parameters, 
m^/M2 and (5, can be obtained to an arbitrary order. 
Technically, the calculations involve integrations of two- 
loop single-scale diagrams of the type of either a muon 
propagator or a massive vacuum diagram. The most 
complicated topology arises when both loop momenta 
are soft, of order m^. This leads to two basic on-shell 
two loop diagrams of the self-energy type which can be 
evaluated using the recent results of Ref. Q . 

The result can be written as 
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and we find 




= — 
48 


= ^ 2 - 


1 ^ ^ 

% dsVsy 

ITtt 


L- 
19 


277 

10368 

Ca 


24^3 3456V3 


128'^^ 288 


IItt^ 

15552 


'M 2 


r 155 2 

( 1 

TT \ 

1296 

V1296 ' 

leVsy 


^L- 


--S2 


+S 


24^3 
1 

^ 72 V 5 J 


36^3 

' A- 


39„ _ 1 . 
64 288^^ 

1 f 


11915 
62208 
3477r2' 
93312 
IItt 


- ^2 -Ca 

384 216^ 


1296 
5377^ "I 


=.52- 


96 ©3 864v^ 


31104 


O 


A2 

M4’ 


( 8 ) 


where L = log(M/m^), (,3 — 1.202057 is the Riemann 
zeta function and S 2 = (f) — 0.260434. 

Substituting a = 1/137.036, Nc = M = 

769 MeV, m^ = 105.66 MeV, = 134.98 MeV, and 
En- = 92.4 MeV into Eq. (0) we obtain the result Q in 
Eq. d). 

The result (i is free from numerical errors, which does 
not mean that it represents the exact contribution of the 
pion pole to < 7 ^ — 2. The dependence of that result on the 
model and form factor parameters has been discussed in 
detail in | 0 |. 

Our result can be checked in several ways, in particular 
to ensure the correct treatment of the Eeynman rules in 
a computer code [|^. Eor example, using the WZW La¬ 
grangian one can evaluate the vacuum polarization con¬ 
tribution to Pp — 2 where the virtual photon splits into 
TT^ and 7 , shown in Eig. 0. The contribution of this dia¬ 
gram to Gp should be positive, since it can be related via 
dispersion relations to the cross section cr(e+e“ ^ 7 r° 7 ) 
&■ Indeed, we find 
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FIG. 2: Vacuum polarization contribution of the neutral pion. 
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Including several more terms in the m^/M and 5 ex- 
pansions, we obtain ^ ~ 3.7 x 10“ , a positive 
contribution. 

It is also useful to compute the pole contribution of a 
(hypothetical) scalar meson “cr” with a mass close to m^. 
By analogy with other diagrams which can be evaluated 
for scalar and pseudoscalar contributions 
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one 
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might expect that the leading ultraviolet logarithm L 
the scalar contribution has equal magnitude but opposite 
sign, compared to the pseudoscalar 7r° |2^ . We have con¬ 
firmed this by an explicit calculation, substituting 
and cr for and 7r° in (||), and we find 
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We see that the leading logarithm is indeed the same as 
in (||) except for the sign. 

Let us also note that the calculation presented in this 
Letter can easily be extended to provide analytic results 
for the pion pole contribution when other models (see e.g. 
iQ) for the pion transition form factors are used. Those 
other form factors can be represented as linear combina¬ 
tions of heavy propagators — and — 
and the technique described in this Letter is applicable. 

Our result confirms the positive sign of the pion pole 
contribution to the muon — 2. We would like to note 
that this effect was found to be positive already several 
years ago in a correct numerical calculation [ p^ . It was 
an apparently simple algebraic mistake which led to the 
sign error in later works . We hope that our analytical 
result will help better understand the dependence of this 
effect on the ultraviolet regulator (M) and thus place the 
theory of the muon g — 2 on firmer footing, as we await 
the release of the new Brookhaven experimental result. 
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